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ABSTRACT. In a previous work by the first author with J. Turi (2), a stochastic variational inequality has been 
introduced to model an elasto-plastic oscillator with noise. A major advantage of the stochastic variational 
inequality is to overcome the need to describe the trajectory by phases (elastic or plastic). This is useful, since 
the sequence of phases cannot be characterized easily. In particular, when a change of regime occurs, there are 
numerous small elastic phases which may appear as an artefact of the Wiener process. However, it remains 
important to have informations on both the elastic and plastic phases. In order to reconcile these contradictory 
issues, we introduce an approximation of stochastic variational inequalities by imposing artificial small jumps 
between phases allowing a clear separation of the elastic and plastic regimes. In this work, we prove that the 
approximate solution converges on any finite time interval, when the size of jumps tends to 0. 

1 . Introduction. The elastic -perfectly-plastic (EPP) oscillator under standard white noise excitation is the 
simplest structural model exhibiting a hysteretic behavior. Moreover, the model is representative of the 
behavior of mechanical structures which vibrate mainly on their first deformation mode. In the context of 
earthquake engineering, relevant applications to piping systems under random vibrations can be accessed 
this way HQ. The main difficulty to study these systems comes from a frequent occurrence of nonlinear 
phases (plastic phases) on small time intervals. A nonlinear phase corresponds to a permanent deformation, 
or in other words to a plastic deformation. A plastic deformation is produced when the stress of the structure 
exceeds an elastic limit. Denoting by x(t) the elasto-plastic displacement, we consider the problem 

x + c x + F(x(s),0 < s < t) = w, (1.1) 



Date: December 21, 201 1. 

"This research in the paper was supported by WCU (World Class University) program through the National Research Foundation 
of Korea funded by the Ministry of Education, Science and Technology (R31 - 20007). 



1 



with initial conditions of displacement and velocity 

x(0) = x , x(0) = y. 

Here Co > is the viscous damping coefficient, k > the stiffness, w is a Wiener process and F({x(s), < 
s < t}) is a nonlinear functional which depends on the entire trajectory {x(s), < s < t} up to time t. 
Denote y(t) := Equation (11.11) written as a stochastic differential equation (SDE) reads 

dy(t) = -(c y(t) + F({x(s),0 < s < t}))dt + dw(t), dx(t) = y(t)dt. (1.2) 

Beyond a given threshold \F({x(s), < s < t})\ = kY for the nonlinear restoring force, the material goes 
through plastic deformation (see e.g. (9l). Introducing A(i), the total plastic yielding accumulated up to 
time t, we can define a new state variable z(t) as z(t) := x(t) — A(t). It follows that in the plastic regime, 
z(t) = 0. From now on, we choose to express the restoring force F({x(s), < s < t}) in (11.21) in terms of 
the new variable z(t) as F({x(s), < s < t}) := kz(t) where \z(t)\ < Y. In other words we consider a 
linear restoring force of the variable z(t) (whose modulus equals Y during the plastic phases). This type of 
force characterizes the elasto-perfectly-plastic behavior. 

In 13, for the previous choice for F, it is shown that (11.21 ) is equivalent to a stochastic variational inequal- 
ity (SVI). In addition, existence and uniqueness of an invariant measure for the solution of SVI have also 
been proven. For a general framework dealing with this class of inequalities we refer the reader to HI and to 
for specific deterministic applications to mechanics. Although SVIs have been already studied in HI to 
represent reflection-diffusion processes in convex sets, no connection with random vibration problems had 
been made so far. From O, the solution (y(t), z(t)) E M 2 of (11.21) satisfies 

y(t) = -(coy(t) + kz(t))+w(t), - y(t))(0 - z{t)) > 0, V|^| < Y, \z(t)\ < Y. (1.3) 

In terms of dynamics of the process (y(t), z(t)), a plastic deformation begins when z(t) reaches and is 
absorbed by Y (resp. —Y) with positive (resp. negative) slope, y(t) > 0. (resp. y(t) < 0) i.e. when 
sign(y(t))z(t) = Y. Then, the plastic behavior ends when the velocity changes sign. At that time, the elastic 
behavior is reactivated. However, around 0, the velocity which is subjected to white noise, changes sign an 
infinite number of times during any small time interval. Often, this leads to a return into plastic behavior in 
a short time duration. This phenomenon is called micro-elastic phasing and has been studied in [81 using 
the numerical method developed in O for the SVI (11.31) . It plays a crucial role on frequency and statistics 
of plastic deformations. Because of this phenomenon, frequency of occurence, statistics (time duration or 
absolute plastic deformation) and the sequence of entry in plastic phase (as well as the sequence of exit) 
are not well defined. In this paper, we consider an EPP oscillator under standard white noise excitation 
subjected to jumps (presented below) to study phase transitions. It has the advantage of separating phases 
clearly, while being an approximation. We prove the convergence of the approximated process towards the 
solution of the stochastic variational inequality (1 1 -3b - 

1.1. Model definition and convergence results. In this subsection, we introduce a stochastic variational 
inequality whose dynamics is "almost" similar to the one of (11.31) except that the second component is 
subjected to jumps of magnitude e > at some random times corresponding to the various exits of the 
plastic phases. 

Precisely, we describe the evolution of the new process (y e (t), z e (t)) by the following procedure; we start 
by defining Tq := and by (y^t), z^t)) the solution of (11.3b . with initial conditions: 

y§(0) = y and 4(0) = z, (y, z) G R x (-Y, Y) := D. 

Then, we define 

t{ := inf{£ > 0, y £ (t) = and \ z e (t)\=Y}. 
For t > t\, let (yl(t), z\(t)) be the solution of (11.31 ) with initial conditions: 

yf(rf) = and z{{rt) = sign^rf)) (Y - e) , 
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again, we define 

r\ := inf{t > t{, yl(t) = and \z\{t)\=Y}. 
In a recurrent manner, knowing r^, and z^(t), we define 

< +1 := inf {t > y«(t) = and |4(i)|=Y}, 
and be the solution of (11.31) with initial conditions: 

Vn+iK+i) = and = sign«« +1 )) (F - e) . 

Now, we define the process (y € (t), z e (t)) on each interval of time [r^, r^ +1 ) as follows: 

y e (t) = -( Co y e (t) + fcz e (t)) + u;(t), (i e (i)-y e (i))(0-z e (i)) > 0, V|0) < F, \z%t)\ < Y (1.4) 
with the following jump-conditions: 

= o, ^«-) = 4-i«), 

and 

= 0, = sign«_ 1 «))(y - 6). 

Remark 1.1. By construction, the process (y e (t) , z e (t)) is cadlag; hence it is regular. In particular, for 
each fixed time T > 0, the number of jumps arise in (0, T], is finite a.s. 

We will prove that the solution (y e (t) , z e (t)) converges to (y(t), z(t)) on any finite time interval, when e 
goes to in the sense described below. 

2. Main results. Our main result is the following theorem. 

Theorem 2.1. Fix T > 0, and consider the processes (y(t),z(t)) and (y e (t) , z e (t)) satisfying d 1 -3b and 
(11.41 ) respectively. Suppose that k > X + (cq) := \ ^ — ^ + c§\J\ + 77zerc, the following convergence 

property holds: 

1 



-E 

e 



sup 

0<t<T 



{| y (t)-^(t)| 2 + fc|z(t)-^)| 2 } 



as e— >0. 



Remark 2.2. Observe that the above condition relating k and cq is purely technical. It will appear clearly 
in the proof of Lemma l25l below. 



2.1. Preliminary results. For (y, z) G D := Kx (— Y, Y), we consider the "elastic" process (y yz (t), z yz (t)): 

~c Q t 1 Cn 1 _32.f+_ \ 

Zyz(t) = e 2 {zcos(w£)H (y H z) sin (a;i)} H — / e 2 i £ s i sin — s))d-u;(s), 

6J 2 uJo 

ft 

2/?/z(*) = z yz(t) + e 2"{— uzsin (ut) + {y + —z) cos (uit)} + / e _ ^~ l _s ' cos (u(t — s))dw(s). 

2 2 Jo 

where, assuming 4 A; > Cq, 

w := — -. 

2 

Remark 2.3. The terminology "elastic" is justified from the observation that (y yz (t), z yz (t)) is actually the 
solution of 

y(t) = -(c y(t) + kz(t))+w(t), z(t) = y(t), (y yg (0), z yz (0)) = (y, z), 

that is the explicit solution of (11.31 ) when the threshold Y = 00 (purely elastic case). Note that the condition 
4k > Cq is needed so that (y(t), z(t)) have real valued solutions. 



Define 

0(y,z) :='m£{t > 0, \z yz (t)\=Y}, (2.5) 

where (y yz (0), z yz (0)) = (y,z). Fort £ [0,T], we set u(y,z,t) := F[9(y,z) > T — t]. This function is 
regular and satisfies the mixed Cauchy-Dirichlet parabolic PDE 

-u t + Au = 0, inD- u{y,Y,t)=0, y > 0; u(y, -V, t) = 0, y < 0; u(y,z,T) = l, (2.6) 

with 

An = ——Uyy + (c y + fcz)itj, - yu z . 

For t < T, the function u(y,z,t) is locally smooth. On the other hand, in the particular case when 
(2/yz(0), %e(0)) := (0, y — e), we consider the probability density function p e of (yo,y-e(*); Zo,K-e(*))- It 
is also known that p e satisfies Chapman-Kolmogorov's equation 

K + ^¥ = 0, p e (y,z,0)=5 0X _ e (y,z), (2.7) 

where A* represents the adjoint operator of A; that is 

A *P £ = -\v yy ~ ((coy + kz)p% + yp\. 

Next, observe that the processes zo Y-e(t) an d yo,y-e(*) are gaussian processes. The key point is to express 
the solution of (12.61 ) through its variational formulation with p e as test function (see proof of Lemma [2~4b . 



The mean, variance and covariance of zo,y- e (t) and yo t Y- e (t) write: 



m e (t) := (Y - e)e~~ (coswt + — smut), a 2 Jt) := / e~ c ° s sm 2 (us)ds, 

2u u l J 

q e (t) := -(Y - e) — e a - sinwt, cr^(i) := / e" c ° s (cos ws - — sina;s) 2 ds, 
w y Jo 2u} 



and 



o-»*(t) := ^e^sin^. 



(2.8) 
(2.9) 

(2.10) 



The density p e then explicitly writes 
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p e (y,z,t) 



cxp 



27ra^)a y (i)(l-p 2 (*)) 1 /2 
_ 2p(*)(y-g e (t))(z-m e (*)) 
a y (t)a z {t) 



2(l-p»(t)) 



(y-g e (i)) 2 , (z-m £ (t)) 



+ 



*2(t) 



(2.11) 



where the correlation coefficient p(t) is defined by (T yz (t)/a y (t)a z (t). Observe that for e = 0. (12.1 II) reduces 
to 



P°(y,z,t) 



exp 



2 TO2 (t) S (i)(l-p 2 (t)) 1 /2 

_ 2p(t)(y-q (t))(z-m°(t)) 
a y {t)a z {t) 



2(1-?®) 



(y-q°(t)) 2 , (z-m°(*)) 



+ 



*3(*) 



with 



n/N c o* , Co. > n/\ t .- c °* ■ 

m (t) := re 2 (coswiH sinwt), g (t) := — Y— e 2 smut. 

*2u to 

From (I2.8I) - (I2.13I) we can easily see that 

m^t) = m°(t) - ef(t), q*(t) = q°(t) + eg(t), 

with 

r I \ — SSL , co. , , N k _ £0l . 

/m := e 2 (coswtH smwt) and g(t) := — e 2 sinwt. 

2w cj 



(2.12) 

(2.13) 
(2.14) 



Plugging (12.141) into (12- lib , we obtain 

p e (y,z,t) 



2n* z (t)<r y (t)(l-f?(t))W 



cxp 



2(1 -p»(t)) 



(y-[g°(*) + gy(*)]) s 



( Z _ [ m o (t) _ e/(t) ] )2 2p(t)(!,- [g°(t) + e^(t)])^- [m°(t)-e/(t)])' 



Oy{t)(J Z (t) 



. (2.15) 



Now, notice that 

(y - [?°(t) + ep(t)]£ + K(t) - e/(t)] ) 2 2p(t)(j, - [q°{t) + e<?(i)] )(z - [m°(t) - ef(t)] ) 



(y - g°(t)) 2 + (z-m°(t)) 2 _ 2p(t)(y- g °(t))(s-m°(t)) + ^ 



2e 



(y-q°(t))g(t) (z-m°(t))f(t) , p(t) 



a y (t)a z (t) 
gf(t) f(t) 2p(t)g(t)f(t) 



Then considering (12.121 ). we have 

P e (y,z,t) =p°{y,z,t)exp 



+ 



((y-q°(t))f(t)-(z-m°(t))g(t)) 



a y (t)a z (t) 

(y - q (t))r(t) - (z - m {t))s(t) 



. (2.16) 



-A(t) + e- 



(1- P 2(t))a y (t)a z (t) 



(2.17) 



where 



A(t) :- 
r(t) := 



g 2 (t) / 2 (t) 2p(t)g(t)f(t) \ 

i-p»(t) v^w ^(o * tf (tMt) y 

s(t)**(t) 



tT-y(t) 



+ P(t)f(t), 



cr z (t) 

We now give a representation of u(0, 1" — e, 0) in terms of the densities p e and p° of the Gaussian processes 
(zo,Y-e{t), yo,Y- e (t)) and (zq : y (t), yo,Y (t)) respectively. The proof is postponed to Section [3] 

Lemma 2.4. Let u be a solution of ( I2.6I ). Then, it satisfies 
u(0,Y - e,0) = 



D 



+ 



L) 



D+ 



[p e {y, z, T) - p°{y, z, T)] dyd 
yu(y,Y,t)p°(y,Y,t) 

r 

yu(y,-Y,t)p (y,-Y,t) 



exD J -l^ A ( t ) + fCjgff) ~ 
expi 2 ^W+ (1 _ p2(0K(tKW 



dydt 



Pi 2 £ AW+ (l-p»(t))a y (t)a,(0 



dydt, 
(2.18) 



with 



Yh(t) := qo ( t )r(t) + (Y-m (t))s(t),h(t):=-g(t)r(t) + (l-f(t))s(t), 
Yl{t) := q (t)r(t)-(Y + m (t))s(t),l(t):=-g(t)r(t)-(l + f(t))s(t), 
D+ := (0, T) x (0, oo) and := (0, T) x (-oo, 0). 



Now consider the terms 

H e -- 



r 



D 



d: 



[p e (y,z,T)-p°(y,z,T)] dydz, 
yu(y,Y,t)p°(y,Y,t) 



cxp 



1 2 A(f) , e(yr(t)-Yh(t)) 

2 U_h (l-f?(t))* v (t)* e (t) 



dydt, 



(2.19) 



yu(y,-Y,t)p°(y,-Y,t) 



1 2, mj e{yr(t)-Yl(t)) 
exp \ --e A(t) + {l _ p2{t))(Jy{t)(Tz{t) 



dydt. 



Next, we study the behavior of these last integrals, with u satisfying ( 12.61 ) so that the previous lemma holds, 
when e is sufficiently small. The proof is also postponed to Section [3] 

Lemma 2.5. Let J € , I e , and H e be the integrals of above. Suppose that 

1 



k > X + (c ) :-- 



3 V 9 6 



Then, 

• lim inf -j = +oo, 

• linig-s.o 4~ is finite, 

• lim^^o is finite. 
Therefore, 

i im "' - r - e -°) 



u(0,-y + e,0) 
+oo ana lim = +oo. 

e->0 



2.2. Proof of Theorem 12.11 We shall use the notation a e n = sign(z e (r^— )). Recall that, for each n > 1, 
the stopping time represents the instant of the n— th jump of the process (y e (t), z e (t)). Hence, for all 
T~n <t < t^ +1 and n > 1, we deduce from d 1 -3b and (11.41 ) that 

y(t) - y e (t) = - [c (y(t) - y e (t)) + k(z(t) - z e (t))} , and 
(z*(t)-y*(t))(z(t)-z*(t))>0, 
(z(t)-y(t))(z e (t)-z(t))>0. 
By using the notation d/dt of derivatives, we obtain 



| (y(t) - y*(t)) = -c (y(t) - y%t)) - k(z(t) - z*(t)), 



(z(t) - z%t)) - (y(t) - y*(t)) (z(t) - z*(t)) < 0, 



(2.20) 
(2.21) 



Multiplying by (y(t) — y e (t)) in (12.201 ) and using the product rule for derivatives, we get from (12.201 ) and 
<HZQ) 



l -^\y{t) - y\t)\ 2 + c,\y{t) - y\t)\ 2 < -k(z(t) - z< (t))(y(t) - y\t)) 



(2.22) 



for all < t < rf l+l and n > 1. Now, integrating (12.221 ) on [r^, t^ +1 ) and noting that y(T^—) = J/(t^), 
y e ( T n~) = y e ( T n) = °. z ( T n~) = z ( T n)> for a11 n > 1, we obtain 



|2/« +1 )| 2 -|y«)| 2 +2c 



\y(t) - tf{t)\ 2 dt+k \z(r* n+1 ) - z^ +1 -)\ 2 -k - z e «)| 2 < 0. 

(2.23) 



But 

k\zK)-Z%T^)\ 2 = fc|(z«)-^«-)) + (^«-)-^«))| 2 

= k\z(^)-z^T^-)\ 2 + ke 2 

+2fce<(z«)-z £ «-)). (2.24) 
Plugging ( 12.241 ) into (12.231 ), and rearranging terms, we obtain 

\yK + i)\ 2 - \yK)\ 2 + k \z(t^ +1 ) - z e « +1 -)| 2 - k |z«) - ^«-)l 2 

f T n+l 

+2c / \y(t) - y £ (t)\ 2 dt < ke 2 + 2fc«z«) - Y). 



We can drop the term 2ke(af l z(T^ l ) — Y) < and get 

\yK + i)\ 2 - \yK)\ 2 + fc H<+i) - ^«+i-)| 2 - & HO - z e «-)\ 2 

+2c / |y(i)-y e (t)| 2 di< A;e 2 . (2.25) 
Observe that, for iV G N*, we can iterate ( 12.251 ) for 1 < n < A" to obtain 

+2c / |y(t)-y e (t)| 2 dt< fciVe 2 . 
Also, recalling that y{r{) = 0, ^(-rf) - z e (r[-)\ 2 = 0, and that /J* |y(i) - y e (i)| 2 dt = 0, we derive: 

1 9 1 9 f T N+l „ „ 

I/C^+i) + fc *(r& +1 ) - z e (r& +1 -) 2 + 2c / |j/(t) - y e {t)\ 2 dt < ke 2 N. (2.26) 

J o 

Denote the total number of jumps of the process (y e (t) , z e (t)) arising in the time interval (0, T) by Nj> := 
maxjv {t^t < T}. Note that T < tjvj+i- Hence, from (12.261 ). we deduce 

sup |y«)| 2 + & sup |z«)-z e «-)| 2 + 2 Co / \y(t) - y*(t)\ 2 dt < ke 2 N^. (2.27) 

l<n<Afj,+l l<n<JV|,+l JO 

Assume first z(0) = Y — e. According to the definition of (12.51 ) set 9 e := 6(0, Y — e) = inf{i > 
0, \z e (t)\ = Y} = ini{t > 0, \z 0)Y -e(t)\ = Y}. It is clear that t{ > 8 e a.s. and then P(r 1 e > 
T) > ¥(6 e > T). Now, let us assume z(0) = -Y + e. It is easy to verify that u(—y, -z, t) = u(y, z, t), 
which gives 

P(0 £ > T) = u(0, Y - e, 0) = u(0, -F + e, 0). 
Thus, by Lemma [231 we have F ^ £ >T ^ — > +oo. Therefore, if the initial condition z(0) associated to (1 1 -4b . 

is a random variable T pi<5y_ e + (1 — pi)5_y+ e independent of the Wiener process w(t), then again 
setting 0f, := inf{t > 0, |*o,r(*)| = F}, 

P(0f > T) 

-> +oo as e -> 0. 

e 

Coming back to equation (12.261 ) and noting that Nj, = ^ n X{t£<t}> we g et 

oo oo 

EA^ = ^E X{t ,<t} = ^X{ r! <T} + ^Ex k <r } . (2-28) 

n=l n=2 

Observe that for all n > 2 and that r £ — r^_ 1 is independent of r £ „ 1 . 



E X{r^<T} = E 



7 



^Kx^^Ex^-^^. (2.29) 



But note that 

From the last inequality and using ( 12.291 ), we deduce 



n-l 



This yields 



EX(r 4S r) < Ex{rE<T}(l - »(0, Y - e,0)) 



E^<Ex W <T, (1 "" ( °- r " e ' 0)) < ^''^ 



Hence, from Lemmas l2.4l and l231 



u(0,Y - e,0) - eu(0,F- e,0) 



e EN? ->0 as 



Thus, as e goes to 0, d2~27l) and (I23TT ) yield 



7 I 



sup |y«)| 2 

l<n<JV|,+l 



+ 2c E / |y(t) - y e (t)| 2 dt + kE 
Jo 



Since the forced jumps have magnitude e, this implies: 



i E 



sup |y(r* 



l<n<N^+l 



fT 

+ 2c E / \y(t) - y e (t)\ 2 dt + £:E 
Jo 



sup 

l<n<JVJ,+i 



sup k«) - 

l<n<iVf,+l 



Also, by (12.221) , we can see that any < t < T^ +1 satisfies 



(2.30) 
(2.31) 

. -»• 0. 
(2.32) 

-> 0. 
(2.33) 



|y(t) - f(t)\ 2 - |y«)| 2 + fc|z(t) - f(t)\* - k \z{r<) - z\r^ < 0. 



f\|2 



This gives 



Hence, 



sup {\y(t) - y%t)\ 2 + fc|z(t) - z*(t)\ 2 } < |y«)| 2 + k |z«) - ^«)| 2 . 

^<*<r* +1 



SUp {|y(t)-^(t)| 2 + A;|z(t)-^(t)| 2 } < SUp 

T{<t<T 



sup 



l<n<JVf+l <<kr< +1 



{l2/(t) - y e (*)l 2 + ^ - ^ e (t)| 2 } j 

< SUp {|y«)| 2 + fc|.«)-^«)| 2 } 



Also, 



sup \\y(t) - ^(i)| 2 + k \z(t) - m\ 2 } < sup {|y«)| 2 + k - ^«)| 2 } . 



0<t<T 

Therefore, (I2.33I ) gives 
1 



-E 



sup 

0<4<T 



\y{t)-y%t)\ 2 + k\z{t)-z%t)\ 2 } 



— > as e— >0. 



3. Proof of the technical lemmas. This section is devoted to the proofs of Lemmas [2.41 and [231 

Proof of Lemma |2~4] From (12.61 ). we have 

= f [ (-ih + Au)p e dydzdt 
Jo Jd 

= 11 (-u t - \u yy + (c y + kz)u y - yu z )p e dydzdt 
Jo Jd 1 

= - p e (y,z,T)dydz + u(0,Y -e,0)+ / / up\dydzdt 
Jd Jo Jd 

^upyydydzdt - J J u((c y + kz)p e ) y dydzdt 



f yu(y,Y,t)p € (y,Y,t)dydt+ [ yu(y, -Y,t)p e (y, -Y,t)dydt 

D~ JD+ 

+ f [ yup e z dydzdt. (3.34) 
Jo Jd 

By using (12.71) and rearranging terms, (13.34b becomes 

u(0,y-e,0)= / p e (y,z,T)dydz+ [ yu(y,Y,t)p e (y,Y,t)dydt- [ yu(y,-Y,t)p e (y,-Y,t)dydt. 
Jd Jd~ Jd+ 

(3.35) 

In addition, p°(y, z, 0) := Sq : y (y, z), and 

0= / p°(y,z,T)dydz+ [ yu(y,Y,t)p°(y,Y,t)dydt - [ yu(y,-Y,t)p (y,-Y,t)dydt. (3.36) 
Jd Jd~ Jd+ 

Using (12.171 ) and substracting (13.361) to (13.351) . we can deduce the result (12- 18b - 

Proof of Lemma |25] First note that, on a neighborhood of t = 0, we have the following expansions: 



CQt 



• f{t) = e~ — (coswt + f^sinwi) = 1 - k\ + f|(cg + 2oj 2 )t 6 + o(t 6 ), 

• g (t) = £e~^ smut = kt(l - ft) + o(t 2 ). 
From (I2~8l)- (I2. 101) . we also have 

. a 2 (t) = f - c ^ + o(t 4 ), = |L ((1 - sft) + o(t)), 

. = ^ (1 - f t + o(t)) , recalling that p(t) = ^f^y ■ 
Equation ( 12.131 ) yields 

f(t) = -Ykt{l-ft) + o{t 2 ), 



m o(t) = y(l-fc4)+o(t 2 )- 



Recalling that r(t) = g( ^ z t) (f) + p(t)f(t) and s(i) := /( * ) ^ ) (t) + and using the previous es- 

timations, we can check that -g(t)r(t) = -*&t + ^c kt 2 + o(t 2 ) and (1 - f(t))s(t) = + 
& (k 2 - | c fc - f ) t 2 + o(t 2 ). Therefore, ~ ^P(c , /c)t 2 where 



P(co,k) := ^ + ^-|, 
9 



Denote X + (cq) := \ \~*3~ + c o y 3 + . Since we have assumed that k > X + (cq), it then follows that 
h'(0) = and h"(0) = ^ p( g°> k ) > 0. We can thus consider a fixed interval (0, t) such that h"(t) > on 



[0, t], hence /i(t) > on (0,t). Also, we have r(t) 



V3, 



(1 — ^P-t + o(t)). Hence, there exists a positive 



constant t such that r(t) > on [0,t]. Let to := min{t,t}. This implies that h(t) > and r(t) > on 
(0, t )- Recall that h(0) = 0. Now write, from d2l9l ) 



with 

n 



t AT r 



yu(y,Y,t)p°{y,Y,t) 



yu(y,Y,t)p°{y,Y,t) 



2- 



2 W (l-p^t))^)^) 



1 



1 



dydt, 
dydt. 



' to AT J -oo 

From the definition of t , h(t) > and r(t) > for < t < to A T. Moreover, y < in if, so we have 



1 2 4W , Miii) 

2 £ W+ (l-p2( t)KWcJz ( t) 



< 0. 



Therefore, the integrand in I\ is a positive function. Now, using the basic inequality exp{— x} — 1 < 
—x exp{— x}, for x > 0, we can write 



> 



t AT r 



1 



yu(y,Y,t)p°{y,Y,t 

e(yr(t) - Yh(t)) 



1 (yr(t)-yfc(t)) 

2 U (1 - pa(t))a y (t)«T.(t) 



xexp ^2 e2 ^ )+ (i-^(t)K(tK(t) 



dydt. 



As A(t) > 0, we get 



> 



t AT r 



yu(y,Y,t)p°(y,Y,t 



(yr(t) - Yh{t)) 
1- p 2 (t))a y (t)a z (t) 
e(yr(t) - Yh(t)) 



xexp ^-2^ + (iV(t)K(tK(t) 



dydt. 



(3.37) 



As the integrand in the right hand side of (13.371 ) is a positive function, Fatou's lemma yields the following 
inequality, 



liminf^ > / / yu(y,Y,t)p°(y,Y,t) 



(yr(t) - Yh(t)) 
(l-pHt))a y (t)a z (t) 



dydt. 



(3.38) 



Note that in (13.381 ) the right hand side may be +00. For I|> since t > to AT, there is no singularity at t = 0. 
Therefore, taking the limit of I|/e, we obtain 



lim inf — 



T r0 



t AT J — oo 



yu(y,y,t)p°(y,y,t)(yr(t)-rM*)) 

(i-p 2 (t)K(tK(t) 



dydt 



(3.39) 



which is finite. Note that 



J 



Kt)| 



(l-p 2 (t))<Ty(t)a z (t) 



yu(y,Y,t)p°(y,Y,t)dy 



dt 
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is finite. Indeed, from the expansion of h(t) we have that locally in time 

Ht) 

<Ty{t)(J Z (t) 

is bounded. Moreover from ( 13.361 ) above 

yu(y,Y,t)p°{y,Y,t)dydt < oo. 



J-oo 



Collecting results we can assert that 



liminf h > f T f iM^Y^Y^ 
e "7o 7-oc (I- p 2 {t))o y (t)a z {t) 



Jo (i-n%(t)^Wv-« / 



The second integral is finite. Now, let us show that the first integral is +oo. We check that 



f 

lim/ y 2 u(y,Y,t)p°(y,Y,t)dy > 0. 
t-+o J_oo 



The function u(y, z, t) is increasing in t. Indeed, from the probabilistic representation we have 

u(y, z,h) = P[0(y, z)>T-h}< P[0(y, z)>T-t 2 }= u(y, z, t 2 ), Vti < t 2 . 
Therefore, we have 

/■0 /-O 

lim / y 2 u(y, Y, 0)p°(y, Y, t)dy < lim / y 2 u(y, Y, t)p°(y, Y, t)dy (3.41) 

Now, 

lty(0-,Y,0) < 0. 

Indeed, Vc > 0, u(-c, Y, 0) > and u(0, Y, 0) = 0. So, 

uj,(0-, Y,0) < 0. 

It cannot be equal to 0, otherwise the derivative exists and u y (0, Y, 0) = 0. But then by minimum properties 
we have u yy (0, Y, 0) > 0, that contradicts 

-u t (o,y,o)--u w (o,Y;o) = o. 

Therefore for y < close to we have 

u{y, Y, 0) ~ ay, a < 0. 
On the interval (—77, 0), we can assume u(y, Y, 0) > |y. So, 



lim f y 2 u{y, Y, 0)p° (y, Y, t)dy > lim f 



° " , n 

y 3 p°(y,Y,t)dy. 



9 

-7} "J-ri* 

From (13.401 ). since ^_ p i^^ ^ a ^ dt = +00, it is sufficient to check the property 

-0 



f 

lim/ y 3 p (y,Y,t)dy < 0. 
3o CO := <?o(0 



Set 

p(i)Y(l - /(t)K(t) 
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then 



P°(y,Y,t) 



Hence, denoting 



2ira y (t)a z (t)(l - p 2 {t))i 



exp 



exp 



L r 



y 



1 



exp 



(y - %(t))' 



we get 

rO 

y 3 p°(y,Y,t)dy 
-r, V2ira z (t) 

In addition, by change of variables, we have 

-gp(t) 



exp 



2(l-p2( t))(J 2 W 

lY 2 (l-/(t)) 2 



(y - %(t)) 2 

2(1 -?Qt))o*(t) J' 



dy, 



du 



(l-p2(t))l/2 CT!/(t) 

Therefore, for t close to we have % (t) 



^^ m) + (1 _ p2 ( )VV,(t)n)3exp(-I^)- 7 =. 



(3.42) 



(3.43) 



Ykt 



and we can check using formula (13.421) and (13.431) that 



lim / y 3 p°(y,Y,t)dy 



J_V3 f° 

2vr 8 



u 3 exp( — u 2 )du. 
2 



Finally, since ^± 
thus have proven 



r(t) _ 3 
(t)cr z (t) ~ W 



(1 + ^f-t + o(t)), the first integral in the right hand side in (13.401 ) is +oo. We 



Next, consider from (12.191 ) the term 



lim — = +oo. 
e->o e 



e 



T roo 




yu(y,-Y,t)p°(y,-Y,t) 



JO 



exp^ 2 ^(t)+ {l _ p2{t))ay{t)(Jz{t) 



dydt. 



From (12.121 ) we have 

P°(y,-Y,t) 



exp 



27ra 2 (t)a y (t)(l-p 2 (t)) 1 /2 
2p(t)(y-q (t))(Y + m (t)) 



2(1 -P 2 {t)) 



(y-q°(t)) 2 (Y + m°(t)f 



*3(t) 



+ 



+ 



a y (t)a z (t) 



Due to the term exp{— 2 (i_p^))^(t) we ^° not nave a singularity because of <r 2 (i). Indeed, Y + m°(i) > 
— exp(— ^j)), and for £ close to we have 



exp 



< 



(y + m°(t)) 2 
2{l - p\t))a 2 z {t) , 

Y + m°(t) 



y(l-exp(-^)) U2(l-pS(W) 



exp 



(y + m°(t)) 2 
2(l-/, 2 (t)K 2 (£) 



< Cexp ( - 



(y + m°(t)) 5 



V 4(l-p 2 (t)V 2 (t), 

where C > is a constant depending on T. From the above equation and recalling the asymptotics a y (t) 
y/i for t close to 0, we deduce that the quantity 

* yu(y, -Y, t)p°(y, -Y, t) L ^I'JjP ,J dydt 




o J o 



L (i-p 2 (£)K(tK(t)- 
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is well defined. In the same way 

— =- I {p e (y,z,T)-p°(y,z,T))dzdy, 

e e JD 

has a well defined limit. Therefore, from (12.181 ). we deduce 

" ( °' y - e ' 0) ^ + oo as £ ^0. (3.44) 
e 

As before, let us assume that z(0) = —Y + e. It is easy to see that u(—y, —z, t) = u(y, z, t). This yields 

F(9 e >T) = u(0, Y — e,t) = u(Q, -Y + e, 0), 
so M (°-~^+ ,: ' ) as W ell. This completes the proof. 
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